Numerical Investigations of Oscillons in 2 Dimensions 
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Oscillons, extremely long-living localized oscillations of a scalar field, are studied in theories 
with quartic and sine-Gordon potentials in two spatial dimensions. We present qualitative results 
concentrating largely on a study in frequency space via Fourier analysis of oscillations. Oscillations 
take place at a fundamental frequency just below the threshold for the production of radiation, 
with exponentially suppressed harmonics. The time evolution of the oscillation frequency points 
indirectly to a life time of at least 10^ oscillations. We study also elliptical perturbations of the 
oscillon, which are shown to decay. We finish by presenting results for boosted and collided oscillons, 
which point to a surprising persistence and soliton-like behaviour. 
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I. INTRODUCTION 

There have been many studies on static non-dissipative 
solutions of classical field theories when the stability of 
the solutions is guaranteed due to a conservation of a 
topological charge (see [lil^) or conserved particle num- 
ber |[3(| and there exist examples where exact solutions 
are known and can be written explicitly in a closed form 
(see e.g. [3, ^Sl). However, energy can be stored in lo- 
calised oscillations of a field for a very long time and in 
this situation there is no obvious conservation law to ex- 
plain the (meta-)stability of the configuration. In 1-1-1 
dimensions 'breather' solutions are known and can be 
written down in a closed form Q. In higher dimensions, 
the long life time of certain oscillations was pointed out 
30 years ago in 0, Q and then re-discovered [9] when the 
dynamics of first order phase transitions and bubble nu- 
cleation was studied. More extended investigations were 
carried out in [10], where conditions, like the size of the 
initial bubble, for long living oscillations were examined. 
In Q oscillating objects were called pulsons. we adopt 
hereafter the definition oscillon according to [lO|. 

There is still not a satisfactory explanation for the long 
lifetime of oscillons. In [ll[ it was suggested that the 
longevity of oscillations could be understand within the 
framework of adiabatic invariance, i.e. there were a con- 
served adiabatic charge explaining the metastability of 
oscillons, called I-balls by the authors [ll|. However, this 
approach assumes negligible gradient energy and non- 
linearity, i.e. quadratic potential. While the adiabatic 
invariance, in spite of its restrictions, provides an ana- 
lytic approach to explain the long, but finite life time, 
some authors have suggested, based on numerical stud- 
ies, that the life time could be arbitrarily long [l^, ll^ . 
and even a Lyapunov exponent was suggested to govern 
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the power law of oscillons life time [13| . 

The possible effects of oscillons on the dynamics of first 
order phase transitions has been studied e.g. in ^1^ (for 
creation of long-living quasilumps from two bubble colli- 
sions see [II])- It was pointed out by Riotto [l^ that at 
the electroweak scale thermal fluctuations are too weak 
to create oscillons and therefore they cannot play a role 
in first order electroweak phase transition. However, a 
very recent study [13] exploited the possibility that oscil- 
lons would speed up the phase transition in the context 
of old infiation. It is natural, that oscillons are not per- 
sistent in a hot heat bath, but they can be long living 
in a weak enough thermal environment (l8| . Principally, 
the necessary conditions of potentials to permit oscillons 
can be fulfilled also in higher dimensional models [l9[, 
but on the other hand it should be noticed that stable 
or extremely long living oscillons (called pseudobreathers 
by the authors) were obtained in sine-Gordon model only 
in two spatial dimensions pol ]. 

In addition to [13] there has been recent interest in os- 
cillons and they have been found in some realistic mod- 
els of great importance. In [IH authors reported that 
with a particular, experimentally relevant, ratio between 
Higgs and W boson masses, there exists an oscillon in an 
SU(2) model that is essentially the bosonic electroweak 
sector of the Standard Model, at zero weak mixing angle. 
An oscillon in the Standard Model could have significant 
consequences e.g. for baryogenesis. Another study [22] 
found oscillon formation after supersymmetric hybrid in- 
flation. The authors worked within the very successful 
D-term inflationary scenario and concluded that for large 
enough, but fully realistic, couplings in the model, oscil- 
lons (called non topological solitons in ) , form in large 
numbers and even dominate the energy density in the 
post-inflationary era (also in (23| energy concentrations, 
often along the domain walls, were reported in simula- 
tions of tachyonic preheating). 

Our study here was motivated by creation of oscil- 
lons from the collapsing domains in theory and in the 
sine-Gordon model, a topic we hope to revisit in another 
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study. Here our investigations concentrate on the prop- 
erties of oscillons when they are created with a Gaussian 
initial ansatz instead of as a consequence of the evolution 
after random initial conditions. However, we try to draw 
the connection by studying the boosted oscillons in the 
end of this paper. Before that we present results from a 
stationary set-up using Fourier analysis to determine the 
time evolution of the oscillation frequency when oscillon 
is created using a spherically symmetric initial ansatz 
and studies of the collapse to spherical symmetry for an 
elliptic ansatz. 

II. THE MODELS AND NUMERICAL SET-UP 

The Lagrangian of a model for a single real scalar 
field (j) in the presence of a potential V is given by 

= (1) 

and the equation of motion thus reads 

'(i>~V^(t) + V'{4,)^Q. (2) 

We study here two potentials with a discrete symmetry: 
degenerate double-well quartic potential 

V{ci,)^\x{^^~v^f (3) 
and sine-Gordon potential 

V{4>) = ^{l-cos{a4^)). (4) 

Scaling out the vacuum expectation value and couplings 
these can be written 

y(</)) = i(^2-lf (5) 

and 

y(0) = -^(l + cos(^0)) (6) 

so that both ^ and ([6]) have minima at </> = ±1 and a 
local maximum at = 0. 

The field equation is evolved on a two-dimensional lat- 
tice with periodic boundary conditions using a leapfrog 
update and a three-point spatial Laplacian accurate to 
0{dx^). The lattice spacing for the data shown is set to 
be dx — 0.25 and the time step dt — 0.05. With that 
choice the fluctuations in total energy on the lattice are 
less than 0.2% over 9 x 10^ iteration steps. We tested 
our method by reducing both dx and dt from the above 
mentioned choice without observing significant difference 
in the quantities of main interest here. Unless otherwise 
stated the simulations for the data shown were carried 
out in 800^ lattices. We choose periodic boundary con- 
ditions rather than absorbing ones because we wish to 
explore the stability of the oscillons to the small pertur- 
bations from any radiation emitted, and because we wish 
to allow them to move without hitting any boundaries. 



III. PROPERTIES 

We start with a Gaussian ansatz 

^(r) =«(l-Aexp(-rVr^)) (7) 

where r is the distance to the center of an oscillon 
r — {x\+X2Y/'^ and the width of the distribution was set 
to be ro ~ 2.9 (in units of (-s/Au)"^), suggested an opti- 
mal choice for a long living oscillon in [lOj] . Earlier studies 
have established the sensitivity both to the spatial size 
ro and the amplitude of the deviation from vacuum. For 
the data shown amplitude A = 1 so that the center of the 
oscillon is located in the local maximum of the potential. 
This choice sets the initial deviation from vacuum to be 
drastic and ensures that we observe the non-linear fea- 
tures of the theory and not a small, linear perturbation 
around the vacuum. 

It should be immediately noted that the Gaussian 
ansatz does not provide the right description of radia- 
tion far from oscillon core. This can be seen by study- 
ing a small, radially symmetric spatial perturbation Lp 
around the vacuum. For quartic potential ([3]) substitu- 
tion 4> = V — ip{r) into ^ leads in lowest order to 

d'^ip(r) 1 dip(r) 9 , x „ 
-^ + — ^ + mVW = 0, (8) 
d'^r r ar 

where = 2Xv^. This Bessel equation has solution 
(f — Ci ■ Ja{mr) + C2 • Yo{mr), where and Y„ are 
the Bessel functions of the first and second kind, respec- 
tively. Both are oscillatory with an amplitude decaying 
asymptotically as r"^''^, thus much more slowly than ([7]). 
Repeating the calculation for sine-Gordon potential dU 
yields the same result with m ~ A^. The parameter m 
defines a mass in the theory as = V"{4>), where (j) is 
at the minimum of the potential. 

Watkins suggested in [l^ the following spherically 
symmetric periodic trial solution (see also |13|) 

00 

(/.(r,t) = ^/„(r)cos(na;i), (9) 

with Lu a free parameter. Studying the coupled ODE's 
for /„ it was noticed that fn fast as n increases, 
so that for n > 5 they are negligible. It was also found 
that solutions existed only for lj < m, going some way to 
explaining the stability: the fundamental mode of oscil- 
lation is not quite fast enough to excite outgoing modes. 
We bear this in mind when discussing our numerical re- 
sults for (/)^ and sine-Gordon potentials. 

A. potential 

The shape of the oscillon profile in cf)^ theory is shown 
in Figure [TJ Deviation from the Gaussian form is clearly 
visible, but not drastic. Figure [5] illustrates the broad- 
ening of the energy density p{x) in the extremum and 
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FIG. 1: Oscillon profile at the extrema (above and below 
the vacuum expectation value) and at crossing the minimum 
of the potential in the 2D 0* theory. Dashed line shows a 
Gaussian form with the same amplitude and the width ro of 
the initial profile. 



contraction to a highly concentrated spike as oscillation 
crosses the minima and kinetic energy dominates the den- 
sity. The value of the field at the centre of the oscillon 
as a function of time is plotted in the Figure [31 An os- 
cillation over one period is not symmetric, as the field 
makes larger excursion and changes more slowly where 
the potential is flatter whereas it changes rapidly in the 
region of steeper potential. Variation of the amplitude 
does not become apparent over this short period of time. 

The total energy and the energy inside shells of radius 
r = 0.5ro, l.Srg and 2.5ro around the centre of the os- 
cillon over the complete simulation time (4.5 x 10^) are 
shown in Figure [H Size r = 0.5ro shows large fluctu- 
ations depending on the phase of the oscillation, which 
can be well understood on the basis of the behaviour of 
p{x) over a period (Figure On the contrary shell 
r = 2.5ro shows relatively thin line and provides a good 
estimator for the energy of an oscillon. This is further 
confirmed by comparing the energy inside shells of ra- 
dius r = 2.5ro and 5ro. They track each other well as 
can be seen in Figure [51 where they are plotted for a 
shorter time interval from the beginning of the simula- 
tion. Figure [5l shows also clearly a period in the begin- 
ning when oscillon radiates energy rapidly. After that the 
rate of energy loss is much weaker and the data suggests 
some evidence for the existence of flat plateaus where en- 
ergy stays at constant value and abrupt drops between 
them. These plateaus could be interpreted as series of 
metastable states (see also [l^l). In the end of the sim- 
ulation more than 70% of the energy is still localised in 
area covering less than 0.5% of the lattice. 

The power spectrum of oscillations was studied by per- 
forming Fourier transforms of (j){t, 0) in consecutive time 
intervals. The length of the interval for the data shown 
is 5000 in natural units, which amounts just over 0.1% of 
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FIG. 2: Energy density p{x) of an oscillon in theory at the 
moment of crossing the minimum of the potential (above) and 
at the moment of maximum excursion (below), corresponding 
the field (f) shown by the lowest solid line in Figure [1] 
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FIG. 3: Field at the centre of the </!>* oscillon 0(t, 0) as a 
function of time t. 
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FIG. 4: Total energy in the lattice and energy inside the 
consecutive shells with radius r = 2.5 ro, 1.5 ro and 0.5 ro 
around the center of the 0* oscillon (from top to bottom) as 
a function of time t. 



FIG. 6: Power spectrum of field at centre of 0'' oscillon [^((.ij)! 
taken in time interval 2.5 ■ 10^ < t < 2.55 • 10^. Solid line is 
a guide to eye for an exponential fit to the amplitudes of the 
first six peaks, exp(— atj/m) with a slope a — 2.03. Vertical 
dotted line shows the radiation frequency u = m. 
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FIG. 5: Detail from Figure |4l showing the total energy in 
the lattice and the energy inside shells with radius r = 5ro 
and 2.5 ro around the center of the </>^ oscillon (from top to 
bottom). 
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FIG. 7: Fit (solid line) to Breit-Wigner form of first peak 
(dots) 0* oscillon taken in time interval 10* < t < 1.25 ■ 10'': 
LOO = 0.8975, r = 0.0020 and K = 0.02. 



the total length covered by the simulation. An interval 
of that length includes approximately 10"^ oscillations. A 
typical example of the power spectrum is shown in Fig- 
ure El There are very distinctive peaks that rise several 
orders of magnitude higher than the backgound between 
them. The first, highest peak just below frequency lo — m 
indicates the oscillation frequency, the other peaks are lo- 
cated at integer multiplies of that. Up to seven peaks can 
be identified in the power spectra. 

There is a significant change in the shape of the peaks 
during the time evolution. In the beginning, when os- 
cillon radiates and the core loses energy reasonably fast, 
peaks are fairly broad, but later on when oscillon 'sta- 
bilises' they become extremely narrow. During the early 



evolution it is possible to make a fit to Breit-Wigner for- 
mula 



cr(w) 



K 



{oj - Lo^f + (r/2)2 ' 



(10) 



where loq is now the peak frequency. A fit to the first 
peak is shown in Figure [T] It cannot be considered very 
accurate, but it gives an order of magnitude estimate of 
the decay width F during the era oscillon radiates its 
energy strongly. Later on the peaks have no width in the 
restricted time interval where the Fourier transformation 
is made (our choice for the length of the interval is fairly 
optimal as longer intervals tend to reveal a shift of the 
peak frequency, not improve the width). 
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FIG. 8; Oscillation frequency corresponding to the location 
of the highest peak as a function of time for 0"* oscillon. The 
precision of the measurement of the frequency is 0.002 at the 
scale of tj/m due to the limited time interval the Fourier trans- 
formation is made. Some data points close to the frequency 
axis are not marked with □ for clarity. Dashed line shows the 
radiation frequency to = m. 



The amplitudes of the peaks in Figure [5] obey clearly 
exponential decay law as a function of the frequency oj. 
The deviations from that at the level shown in the figure 
are most likely once again due to the limited accuracy of 
discrete Fourier transformation. The data also points out 
a low frequency peak (just next to |<?!)(a;,0)| axis), with 
frequency approximately 0.06. There is a slight variation 
of the amplitude, a beat, with the corresponding period 
around 100 in time units (though not present in Figure[31), 
which we believe is the cause of the structure. 

Not only the shape of the peaks but also their loca- 
tion changes during the time evolution. There exists 
a strong correlation between the energy in the oscillon 
(shell radius r — 2.5ro in FigureS]) and the oscillation fre- 
quency, i.e. location of the highest peak, which is shown 
in Figure [8] as a function of time. During the early pe- 
riod when the oscillon radiates strongly, the frequency 
increases rapidly, but then the growth slows down dras- 
tically as the rate of energy loss becomes tiny. The time 
evolution of the oscillation frequency is the key to pre- 
dict the life time of an oscillon: as long as the frequency 
stays below radiation frequency, oscillon cannot directly 
radiate all of its energy and disappear. As the rate of en- 
ergy loss decreases, the increase in frequency slows down. 
Turning it other way round this is in agreement with 
Watkins [13 who reported the radiation rate decrease as 
Lj/m increases, with a minimum at w/m ~ 0.97. 

Though we have not evolved oscillons in our simula- 
tion longer than 4.5 x 10^ time units, we can give bounds 
on the life time of oscillon on the basis of the evolution 
of the oscillation frequency. Even a linear fit to second 
half of the points in Figure [5] suggests that the radiation 
frequency is not reached before a time of a few times 10^, 
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FIG. 9: Oscillon profile at the extrema (above and below the 
vacuum expectation value) and at crossing the minimum of 
the potential in the 2D sine-Gordon theory. Dashed line shows 
a Gaussian form with the same amplitude and the width ro 
of the initial profile. 



or about 10'' oscillations. Similar lifetimes were reported 
in (2^ where lattices with a technique of adiabatic damp- 
ing were exploited. However, the lifetime can be much 
longer as the slope of the line in Figure [5] is clearly flat- 
tening out as well as the strongly decreasing radiation 
rate as a function of frequency reported in jl2l |. 

Setting the original amplitude A too high (e.g. A = 10 
to the direction where the potential is steeper) will not 
lead to an oscillon, the initial concentration of energy 
does not stay localised but spreads rapidly. Starting with 
a configuration that has more energy than the one shown 
in Figure [1] can still evolve to an oscillon but with a mod- 
ulated amplitude. This effect becomes most apparent in 
the case of sine-Gordon potential. 

B. Sine-Gordon potential 

Persistent oscillations in radially symmetric sine- 
Gordon equation were studied in [13] and extreme sensi- 
tivity both to the form and the amplitude of the initial 
deviation from the vacuum were observed. A complicated 
initial profile was used to obtain a minimally radiating 
pseudobreather. We have not been able to create a 'min- 
imal' oscillon with a Gaussian initial ansatz - though the 
profile of an oscillon in sine-Gordon potential is also fairly 
close to Gaussian (Figure [51) oscillations have modulated 
amplitude fFigure fTU)) the period of which depends on the 
initially set amplitude A. Figure [TT] shows total energy in 
the lattice and energy inside shells of several radii around 
the centre of the oscillon in a simulation that spans the 
evolution up to time 10^. There is no apparent evidence 
for energy plateaus in Figure [TTl 

The same choice of the interval for Fourier transfor- 
mation was used as in 0"* theory that now corresponds 
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FIG. 10: Field at the centre ol the sine-Gordon oscillon 
(0(t,O)) as a function of time t. 
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FIG. 11: Total energy in the lattice and energy inside the 
consecutive shells with radius r = 5.0 ro, 2.0 ro and 1.0 ro 
around the center of the sine-Gordon oscillon (from top to 
bottom) as a function of time t. 



0.5% of the total time of the simulation. The study in 
frequency space shows similar basic features as in <j>'^ the- 
ory with some significant differences. The sine- Gordon 
potential ^ is symmetric around minima in contrast 
to quartic potential ([3]). Due to this symmetry there 
are no even harmonic peaks in the power spectrum of 
sine- Gordon oscillon as can be seen in Figure [T^] and at 
very best only first four peaks were visible in the power 
spectrum. As the amplitudes of the peaks decrease ex- 
ponentially also here, but only half of them are present 
compared to Figure [5] of cj)'^ theory, this suggests that the 
ansatz ^ would be a particularly good approximation 
for an oscillon in sine-Gordon potential, or potentially 
in any other symmetric potential that allows oscillons. 
There is no sign of low frequency structure in Figure [T^] 
in contrast to Figure [S] 



FIG. 12: Power spectrum of field at centre of sine-Gordon 
oscillon taken in time interval 1.5 ■ 10^ < t < 1.55 • 

10^. Solid line is guide to eye for an exponential fit to peak 
amplitudes, exp(— acj/m) with a slope a — 1.84. Vertical 
dotted line shows the radiation frequency to = m. 
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FIG. 13: Oscillation frequency as a function of time for sine- 
Gordon oscillon. The precision of the measurement of the 
frequency is 0.003 at the scale oi uj/m. Dashed line shows the 
radiation frequency w = m. 



The time evolution of the frequency of the first peak 
in power spectrum of sine- Gordon oscillon is shown in 
Figure [T3] from a simulation that was evolved 10^ time 
units. The oscillation frequency starts initially closer to 
the lowest radiation frequency compared to the quartic 
potential, but there is no drastic increase either, which 
may be due to the weaker coupling to radiative modes 
as even multiples of the oscillation frequency are absent. 
Also here the increase in the oscillation frequency slows 
down in the course of time and even linear extrapolation 
yields a life-time estimate, i.e. intersection with the ra- 
diation frequency, around 10^ time units, but it could be 
much larger, especially as the study of the power spec- 
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FIG. 14: EUipticity, |Oii|, as a function of time in the 2D 
theory. Straight, solid line is guide to eye of an exponential 
fit, exp{—bt) with a slope b = 0.009. 



trum suggests weaker radiative modes. 

IV. ELLIPTICAL OSCILLONS 

The relative ease to create an oscillon together with 
the energy shells pointing to the existence of successive 
plateaus and thus a series of metastable states gives rise 
to the question if the field configuration of an oscillon is 
unique. A very modest approach to the uniqueness is a 
study of an initial configuration which is not spherically 
symmetric. As our numerical method was not restricted 
to spherical symmetry a Gaussian ansatz with different 
width in xi and X2 direction was used. Oscillons emerg- 
in g fro m the collapse of asymmetric bubbles was studied 
in [2^ using effective radius as a measure. We study here 
the time evolution of the ellipticity via the quadrupole 
moment gij of the energy density x) 

Jd^xp{t,x) ■ ^ ' 

Ellipticity is given by the eigenvalues of the traceless ma- 
trix 

'di] ^ Qij -^SijQ.j . (12) 

Because the major and minor axis are along xi and X2 the 
off-diagonal entries 812, 621 vanish (a good check for the 
numerical accuracy of the method: |0i2|, |62i| < 10"^'"') 
and then ellipticity is given directly by the diagonal el- 
ements of ([T^ . We measured Qij in a square of length 
20 in physical units located in the center of the lattice. 
Lattice size was set to be 4800^ and thus no boundary 
effects will be important until far after time = 600. Fig- 
ure [U shows |9ii| as a function of time when initially 
the ratio of the major to the minor axes was set to be 




FIG. 15: Contours of 0(x) for the initial elliptic Gaussian 
profile and after two oscillations (time = 10.6) in the 2D 0'* 
theory. 
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FIG. 16: Total energy and energy inside shells of radius r = 
2.5 ro and r = 1.0 ro around the center of oscillon as a 
function of time (from top to bottom). Here ro = 3.9 is the 
initial width of the major axis. 
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FIG. 17; EUipticity, |Oiil, as a function of time in the 2D 
sine-Gordon model. Straight, solid line is guide to eye of a 
power law fit, with a slope S = 2.87. 
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FIG. 18: Total energy and energy inside shells of radius r = 
2.5 ro and r = 1.0 ro around the center of sine-Gordon oscillon 
as a function of time (from top to bottom). Here ro = 3.9 is 
the initial width of the major axis. 



3 : 2 (with the width in the xi direction being ro = 3.9). 
The oscillon approaches rapidly a spherical profile and 
ellipticity disappears in exponential phases. The orien- 
tation of the major axis chances in a way that is not ob- 
viously periodic (zero crossings of Gn appearing in the 
Figure [Ml as almost vertical lines, which is demonstrated 
much more illustratively in Figure llSp . Figure shows 
total energy and energy inside two shells around the os- 
cillon core. Though spherically symmetric shells are not 
an entirely adequate method to measure energy of an el- 
liptic configuration, they give a reasonable estimate and 
show that while ellipticity disappears exponentially de- 
creasing four orders of magnitude, the energy inside the 
oscillon core has decreased by just over 10% indicating a 
quick collapse to a spherical shape. 



Figure fTTI shows the ellipticity and [TSl total energy and 
energy inside shells in a simulation with the same set-up 
but for sine-Gordon potential. Qualitatively features are 
similar, ellipticity decays rapidly and oscillon collapses 
into spherical profile. Here the orientation of major axis 
changes frequently in almost periodic way. Moreover, the 
decay of ellipticity |8ii| (in Figure [T7)l matches better to 
a power law than an exponential fit. It is very tempting 
to try to interpret this result on the basis of the obtained 
differences of oscillons in frequency space between quartic 
and sine-Gordon potentials. As both even and odd mul- 
tiples of the oscillation frequncy are present in case of 
the quartic potential, the first radiative frequency, twice 
the basic frequency, has far greater amplitude than in 
sine-Gordon potential where the first radiative mode has 
stronger suppression its frequency being three times the 
basic oscillation frequency. Therefore it seems plausible 
that oscillon in quartic potential can radiate its asymme- 
try exponentially, while initial deviation from spherical 
symmetry decays only with a power law in sine-Gordon 
potential. Naturally, the study above does not exclude 
the possibility of different spherically symmetric forms 
for an oscillon. 



V. COLLIDING OSCILLONS 

If created in Nature via some mechanism like collapse 
of domains, oscillons would have initially translational 
momentum. In order to draw the connection between os- 
cillons given birth in simulations with random initial con- 
ditions and stationary oscillons created with e.g. Gaus- 
sian ansatz, moving oscillons were prepared. Practically 
the oscillon that starts with a Gaussian initial profile is 
first allowed to evolve few oscillations and then the con- 
figuration is Lorentz boosted. For all the shown data the 
velocity of moving oscillons is set to be 0.5 (we point out 
that oscillons originating from collapsing domains seem 
to generally have higher velocities). 

Fast moving oscillons become some sort of waves: they 
have typically at least two 'phases', part of the wave front 
lies above, part below the minimum of the potential. 
There is a precise correspondance between boosted os- 
cillons and the objects travelling on wave fronts originat- 
ing from collapsing domains. We have also checked that 
boosted oscillon is persistent like a stationary one and 
can circulate on a lattice (with periodic boundary condi- 
tions) travelling long distances without demise. Here it is 
appealing to speculate that there is a connection to freak 
waves in the context of pattern formation in oceonogra- 
phy (for a study of freak wave formation from a nonlinear 
Schrodinger equation see [26i|). 

We finish by reporting the results of studies where os- 
cillons were made to collide. There is twofold reason- 
ing for this study. We want to examine what kind of 
interactions oscillons have and thus try to understand 
their nature from that point of view. Secondly, if created 
in large numbers in phase transitions, collisions are in- 
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FIG. 19: Sequence of snapshots at times t — 0, 36.25, 80 of 
an on-axis collision of two oscillons, an immobile (originally 
located in the middle) and moving (from right to left), in the 
2D theory. 



FIG. 20: Sequence of snapshots at times t = 0, 36.75, 100 of 
an on-axis collision of two oscillons in the 2D theory. The 
phase at collision set 1/4 of the period different compared to 
Figure [H 



evitable and the behaviour of oscillons in those has thus 
impact on the subsequent evolution of the system. The 
general rule is that a collision does not destroy an os- 
cillon: we have never witnessed that the first collision 
between oscillons would have lead to a demise of either 
of them though it is noticeable that oscillons disturbed 
this way have also disappered in the second collision that 
occured on a lattice with periodic boundary conditions. 

Figure [19] shows snapshots from a collision between im- 
mobile and moving (velocity 0.5) oscillons. After the col- 
lision there is one fast moving oscillon continuing to the 
direction of original translational momentum and fairly 
stationary oscillon whose center, however, has moved 
slightly from the original location. This example seems 
thus almost interactionless. However, amount of inter- 
action and momentum transfer depends on the phases 



of oscillons at the moment of collision and can result 
to two moving oscillons. This can be examined as the 
set-up for collision has been prepared by creating the im- 
mobile oscillon at arbitrary time on the lattice, we can 
control the phase of stationary oscillon via the moment 
of its creation. The effect of the phase is demonstrated 
in Figure [201 where the set-up is otherwise the same as 
in Figure [TOl but the initially stationary oscillon is time 
1.3 behind the previous one corresponding approximately 
1/4 of the period. Here both the oscillons move left after 
the collision. 

Figure [51] shows snapshots of a coUison between two 
boosted oscillons with the same velocity (0.5) travelling 
into opposite directions. As can be seen oscillons pass 
through each other. We have also changed the align- 
ment of oscillons so that the collision does not occur head 
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FIG. 22: Sequence of snapshots at times t = 0, 22, 49 of an 
off-axis collision of two oscillons in the 2D 0^ theory. 



on. Figure [22l shows snapshots were the deviation in the 
aUgnment between the centers of oscillons is 5.0 in phys- 
ical units. The set-up leads to an attractive scattering, 
oscillons do not behave as classical particle like objects, 
but their paths bend towards each other the angle be- 
tween initial and final velocity being approximately 20°. 

FIG. 21: Sequence of snapshots at times t = 0, 29.5, 
30.5, 31.5, 60.5 of an on-axis collision of two oscillons in the 
2D <?i* theory. 
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VI. CONCLUSIONS 

We have studied oscillons, extremely long living oscil- 
lations of a scalar field, numerically in two dimensions 
for and sine-Gordon potentials. We evolved station- 
ary oscillon for 10^ time units or longer (which is three 
orders magnitude larger than the reported life time for 
I-balls, lO'^m"^, in [ll|) and examined the power spec- 
trum of oscillations. Study in the frequency space not 
only suggested much longer life time on the basis of the 
time evolution of the basic oscillation frequency than we 
could directly probe, but also points out the validity of 
the separable ansatz Though we have not evolved 
the oscillon in sine-Gordon potential as long, the analy- 
sis of the power spectrum suggests that it might be even 
more persistent than oscillons in cfA theory in spite of the 
modulated amplitude. Studies of Lorentz boosted oscil- 
lons and collisions between them show also persistence 
of these objects and point out the behaviour that is very 
similar to solitons. 

As our study was carried out on lattices with periodic 
boundary conditions, the radiation oscillons emit comes 
to some extent back to the centre of the lattice where 
oscillon is located. This can be viewed as a small per- 
turbation and long survival of oscillons in the set-up as 
their persistence to radiation. On the other hand incom- 
ing radiation can also work otherwise pumping en ergy 
back into oscillon and thus extending the life time [1^ . 
Even in the latter case our study has relevance in situa- 
tion where oscillon can absorb energy from enviroment, 
like a weak heat bath or a laboratory experiment in a 
closed system. 

There are several unanswered question concerning os- 
cillons and their longevity. The mechanism that makes 
some potentials able to confine energy over millions of os- 
cillations should be understood together with the relation 
to the oscillation frequency. On the basis of the spreading 
of the energy density during the oscillation (see Figure [2]) 
one would naively assume that energy could escape in 
the form of radiation relatively quickly. The oscillation 
frequency also seem to be fined tuned just below the ra- 
diation frequency. This has been reported also for the 
SU(2) model oscillon [21] where the difference between 
oscillation and radiative frequencies was observed to be 
per mil level. Finally, even the level and role of nonlin- 



earity in the nature of oscillons is not well understood. It 
has been noticed, both in the study of freak wave forma- 
tion from a nonlinear Schrodinger equation (26j as well 
as oscillon (called axiton by the authors) formation in 
the axion field [27,], that nonlinear effects are crucial in 
the generation of high density contrasts. On the other 
hand, in TT] the authors emphasise the importance of 
the potential being dominated by the quadratic term as 
a condition for I-balls to exist as well as in [l2| weak- 
ness of nonlinear effects was considered important for the 
longevity of oscillons. 

Almost stable, non-dissipative solutions of field theo- 
ries oscillons have interest itself, but also their existence 
can lead to significant consequences that are worth stud- 
ies of their own. In particular, if formed, oscillons could 
play a role in the early universe and phase transitions 
that occured there especially if they survive long enough 
in realistic conditions, as for instance their behaviour 
in collisions suggests. Here e.g. the oscillon in sine- 
Gordon potential may have relevance in the QCD phase 
transition as the axion potential has sine-Gordon form. 
Kolb and Tkachev have studied the formation of solitons 
in [27[. While finishing this paper we became aware of 
the study [1^ where similar life times of oscillons as we 
obtained here were reported, but in an expanding back- 
ground in one dimension with a potential including also 
(f)^ -term. If this longevity persists in full three dimen- 
sional models with Hubble expansion, oscillons may have 
interesting cosmological implications. 

Note added: while this paper was being considered for 
publication the study [291 ] appeared examining oscillons 
in 3 dimensions. 



Acknowledgments 

The authors thank Nick Manton, Bernard Piette and 
Anders Tranberg for discussions and Paul Saffin for 
pointing out the reference [T^ . P. S. was supported 
by Marie Curie Fellowship of the European Community 
Program HUMAN POTENTIAL under contract HPMT- 
CT-2000-00096 and by the Netherlands Organization for 
Scientific Research (N.W.O.) under the VICI programme 
and acknowledges the hospitality at the University of 
Sussex. 



]1] G. 't Hooft, Nucl. Phys. B79, 276 (1974). 

[2] A. M. Polyakov, JETP Lett. 20, 194 (1974). 

J3] S. R. Coleman, Nucl. Phys. B262, 263 (1985). 

[4] M. K. Prasad and C. M. Sommerfield, Phys. Rev. Lett. 

35, 760 (1975). 
[5] S. Theodorakis, Phys. Rev. D61, 047701 (2000). 
]6] R. Rajaraman, Solitons and Instantons (North-Holland, 

1982). 

[7] L L. Bogolyubsky and V. G. Makhankov, JETP Lett. 24, 
12 (1976). 



[8] I. L. Bogolyubsky and V. G. Makhankov, Pisma Zh. 

Eksp. Teor. Fiz. 25, 120 (1977). 
[9] M. Gleiser, Phys. Rev. D49, 2978 (1994), hep- 

ph/9308279. 

[10] E. J. Copeland, M. Gleiser, and H. R. MuUer, Phys. Rev. 

D52, 1920 (1995), hep-ph/9503217. 
[11] S. Kasuya, M. Kawasaki, and F. Takahashi, Phys. Lett. 

B559, 99 (2003), hep-ph/0209358. 
[12] R. Watkins, DART-HEP-96/03, unpublished (1996). 
[13] E. P. Honda and M. W. Choptuik, Phys. Rev. D65, 



12 



084037 (2002), hcp-ph/0110065. 
[14] M. Glcisor and R. C. Howell (2003), cond-mat/0310157. 
[15] 1. Dymnikova, L. Koziel, M. Khlopov, and S. Rubin, 

Grav. Cosmol. 6, 311 (2000), hep-tli/0010120. 
[16] A. Riotto, Phys. Lett. B365, 64 (1996), hep-ph/9507201. 
[17] M. Gleiser (2006), hep-th/0602187. 
[18] M. Gleiser and R. M. Haas, Phys. Rev. D54, 1626 (1996), 

hop-ph/9602282. 
[19] M. Gleiser, Phys. Lett. B600, 126 (2004), hep- 

th/0408221. 

[20] B. Piette and W. J. Zakrzewski, Nonlinearity 11, 1103 
(1998). 

[21] E. Farhi, N. Graham, V. Khemani, R. Markov, and 
R. Rosales, Phys. Rev. D72, 101701 (2005), hep- 
th/0505273. 

[22] M. Broadhead and J. McDonald, Phys. Rev. D72, 



043519 (2005), hep-ph/0503081. 
[23] E. J. Copeland, S. Pascoli, and A. Rajautie, Phys. Rev. 

D65, 103517 (2002), hcp-ph/0202031. 
[24] M. Gleiser and A. Sornborger, Phys. Rev. E62, 1368 

(2000), patt-sol/9909002. 
[25] A. B. Adib, M. Gleiser, and C. A. S. Almeida, Phys. Rev. 

D66, 085011 (2002), hep-th/0203072. 
[26] M. Onorato, A. R. Osborne, M. Serio, and S. Bertone, 

Phys. Rev. Lett. 86, 5831 (2001), cond-mat/0104055. 
[27] E. W. Kolb and I. I. Tkachev, Phys. Rev. D49, 5040 

(1994), astro-ph/9311037. 
[28] N. Graham and N. Stamatopoulos, Phys. Lett. B639, 

541 (2006), hep-th/0604134. 
[29] G. Fodor, P. Forgacs, P. Grandclement, and I. Racz 

(2006), hep-th/0609023. 



